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Abstract 

Ah initio potential energy and dipole moment surfaces are computed for sulfur trioxide (SO3) at 
the CCSD(T)-F12b level of theory and appropriate triple-zeta basis sets. These surfaces are fitted 
and used, with a slight correction for the equilibrium S-0 distance, to compute pure rotational and 
rotation-vibraton spectra of ^^S^^Oa using the variational nuclear motion program TROVE. The 
calculations considered transitions in the region 0-4000 cm~^ with rotational states up to J = 85. 
The resulting line list of 349 348 513 transitions is appropriate for modelling room temperature 
^^S^^Oa spectra. Good agreement is found with the observed infrared absorption spectra and the 
calculations are used to place the measured relative intensities on an absolute scale. A list of 10 878 
experimental transitions is provided in a form suitable for inclusion in standard atmospheric and 
planetary spectroscopic databases. 
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I. INTRODUCTION 



Sulfur trioxide (SO3) is a stable, planar, symmetric molecule whose electronic ground 
state is a closed shell. On earth it is a pollutant produced from smoke-stacks and other 
industrial exhausts [1]. In the atmosphere SO3 forms sulfuric acid with its association with 
acid rain, and inside combustion systems it is corrosive. In both cases the reactivity of SO3 
makes it production hard to monitor j^. SO3 is produced naturally on earth as part of 
volcanic emissions [3] . It is also thought to be a significant constituent of the atmosphere of 
Venus Q]. 

The infrared vibration-rotation spectrum of ^^S^^03 (henceforth referred to as SO3) has 
been extensively investigated in a series of papers by Maki and co-workers jJInl. Its "forbid- 
den" rotational spectrum, for which centrifugal distortions can induce transitions, has been 
investigated using microwave Fourier-transform spectroscopy However, although these 
works provide an extensive list of measured line frequencies, none of them report absolute 
transition intensities. 

There has been limited theoretical work on SO3. Dorney et al {isl reported force 
constants, while Martin [14] computed an ab initio potential energy surface (PES) using 
coupled cluster methods and reported theoretical estimates for the band origins of the low- 
lying vibrational states. Again, neither of these works considered transition intensities. 

The lack of any absolute transition intensities for SO3 places severe limitations on the use 
of its infrared spectrum for remote sensing; applications or inclusion of this data in standard 

n n 

atmospheric and planetary spectroscopic databases [15|, Hq. In this work we report the 
calculation of a new ab initio PES and associated dipole moment surfaces (DMS) for SO3. 
These are used to not only produce theoretical spectra for ^^S^^03, but also to place the 
relative intensity measurements of its infrared spectrum on an absolute scale. 

II. THE POTENTIAL ENERGY AND DIPOLE MOMENT SURFACES 

The ab initio PES was computed using the recently-proposed explicitly correlated F12 
singles and doubles coupled cluster method including a perturbational estimate of connected 
triple excitations, CCSD(T)-F12b [ill, conjunction with the corresponding F12-optimized 
correlation consistent basis sets, namely the valence correlation-consistent functions aug-cc- 
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pVTZ-F12 and aug-cc-pV(T+d)Z-F12 for oxygen and sulfur, respectively |l8| . We also 
utilized the OptRI Q cc-pV5Z/JKFIT and aug-cc-pwCV5Z/MP2FIT lE] auxiliary 
basis sets for evaluating the many-electron integrals, the exchange and Fock operators, and 
the remaining electron repulsion integrals, respectively. The value of the geminal Slater 
exponent /3 was chosen as 1.2. In ab initio calculations Molpro.2010 [22] was employed. We 
used a large grid of 13 000 points with the stretching modes ranging from 1.1 to 2.0 A, 
interbond angles a between 70 and 123°, and the inversion angle p between and 50°. This 
grid was sufficient to cover energies up to 40 000 cm~^ above the minimum. 
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This PES was then represented in an analytical form given by the expansion 

V (6 ,6,6, 6a , 66 ; sin p) = Ve + Vo (sin p) + ^ Fj (sin p) 

j 

+ XI ^■'■'^ ^^^^ ^) 6 + XI ^^''^ ^^^^ ^' 

j<k 3<k<l 

j<k<l<'m 



in the coordinates ^k'- 



where 



and 



= 1 - exp(-a(rfc - re)), fc = l,2,3, (2) 

^Aa = (2q^i - ^2 - Ol'i) , (3) 

^46 = ^ (^2 - ^^3) , (4) 

2 

sinp = — sin[(ai + 0^2 + <^3)/6], (5) 
V3 



1/0 (sin p) = ^ /(S'^ (sinpe - sinp)^ (6) 



Fjk.S^mp) = ^ /jj (sinpe - sinp)^ (7) 
The same form has been used to represent the PESs of NH3, PH3, SbHa, NH3+, and BiHa 



27| . The potential parameters /-^ were obtained through a least squares fit to the ah 



initio points with an root-mean-squares (rms) error of 0.067 cm ^ . Weight factors were set 



using the expression suggested by Partridge and Schwenke [28|]: 



tanh [-0.0005 cm x {V^ - 16 000 cm-^)] + 1.002 002 002 
2.002002002 cm X Vp 
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where V^'^^ = max(16 000 cm~-^, 1/^), and Vi is the ab initio energy at the ith geometry (in 
cm~^), measured relative to the equihbrium energy. The ab initio energy Vi is weighted 
by the factor Wi in the PES fitting; these weight factors favor the energies below 16 000 
cm~^. The ab initio equilibrium geometry Tq as obtained from the fitting is 1.42039 A. This 
can be compared with an experimentally derived value 1.41732 A [6]. We found that the 
experimental rotational energies of SO3 [8] are best described by the latter value. Therefore 
we decided to use this value of in all our calculations in place of the ab initio value. It is 
known however that replacing the equilibrium structure may cause undesirably large changes 
to the ro- vibrational energies. In order to minimize this effect the following procedure was 
employed. 

The idea is to expand the ab initio PES around some reference value Tref and than 
remove the linear terms in the PES expansion and replace the expansion center with re^^^^ . 
In doing this we expect the change of the shape in the PES to partly compensate the effect 
from replacing re^^^ with re^^^^ on the vibrational energies. We choose the reference center 
r^^^^^ = 1.42039 + Are, where Ar^ = re^^^ — re^^^^ = 0.00307, i.e. on the opposite side from 
re^^^^ and at the same distance from ri^^'^K 

This is the only adjustment to the shape of ab initio PES utilized. Table [J compares 
rotational levels obtained by Maki et al. [8] and those computed using the equilibrium- 
adjusted PES used in this work. The agreement is very good. 

The DMS were calculated using the same level of theory as the PES and on the same 
grid of 13 000 points. The ab initio values were then expressed analytically using the 
symmetrized molecular bond (SMB) representation described in detail in Ref. [29|]. The 
resulting dipole moment parameters obtained through a least squares fit reproduce the ab 
initio data with an rms error of 0.00013 D. In these fittings the same factors defined by 
Eq. ([8]) were used to weight the geometries according to the corresponding energies. To our 
knowledge there are no experimental or ab initio dipole moment data in the literature that we 
could use to validate our DMS against. However, our experience of working dipole moments 
for different systems (29l-l32[| shows that ab initio intensities in most cases are competitive 
with experimental measurements. The quality of the relative intensities calculated using our 
DMS is discussed below. 

Both the potential energy and dipole moments functions used in the present work are 
given as supplementary data. 
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TABLE I: Theoretical rotational term values for ^^S^^Oa ground vibrational state (cm ^) compared 
with experiment [s]. 



J 


K 


Obs. 


TROVE 


2 





2.0912 


2.0916 


3 


3 


2.6115 


2.6119 


4 


3 


5.3998 


5.4006 


4 





6.9707 


6.9718 


5 


3 


8.8852 


8.8864 


6 


6 


8.3548 


8.3559 


6 


3 


13.0675 


13.0694 


7 


6 


13.2342 


13.2360 


7 


3 


17.9467 


17.9493 


8 


6 


18.8106 


18.8132 


8 


3 


23.5228 


23.5263 


8 





25.0935 


25.0972 


9 


9 


17.2297 


17.2319 


9 


6 


25.0838 


25.0874 


9 


3 


29.7958 


29.8002 


10 


10 


17.2297 


17.2319 


10 


9 


24.2002 


24.2035 


10 


6 


32.0539 


32.0584 


10 


3 


36.7655 


36.7709 


10 





38.3360 


38.3417 


20 


18 


89.8252 


89.8372 


20 


15 


107.0973 


107.1122 


20 


12 


121.2253 


121.2425 


20 


9 


132.2114 


132.2304 


20 


6 


140.0574 


140.0777 


20 


3 


144.7645 


144.7856 


20 





146.3334 


146.3548 


80 


78 


1195.6589 


1195.8085 



III. RO- VIBRATIONAL CALCULATIONS 

A. Basis set convergence and Hamiltonian optimisation 



Ro-vibrational calculations were performed with the program TROVE j33| adapted to 
work in the V^i^{M) molecular permutation-inversion group appropriate for SO3. In order 
to achieve results of high accuracy as well as minimising the requirement for computational 
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resources, it is necessary to optimise the size of the Hamiltonian matrix. This involves 
prehminary truncation of the basis set, as weU as hmiting the order of both the kinetic and 
potential components of the Hamiltonian expansion. TROVE employs a polyad number 
truncation which controls the size of the basis set. For SO3 the polyad number is given by 

P = 2(ni + 712 + ns) + n4 + ns + y , (9) 

where rii are the quanta associated with ID basis functions, 0^, whose product gives our 
vibrational basis set [33]. Each of these basis functions is associated with an internal co- 
ordinate and only functions for which P < Pmax are included in the primitive basis set. 
Initial tests were carried out to measure the degree of convergence using different values for 
^max, and the orders of the kinetic and potential energy expansions. In this work we use 
a kinetic energy expansion of order 4, and a potential energy expansion of order 8; using a 
kinetic energy expansion order of 6 requires a more expensive calculation where convergence 
is already observed to within 0.001 cm~^ when expanding to fourth order. In the present 
study, we find that the convergence is more sensitive to Pmax, and we obtain convergence to 
within 0.1 cm~^ when Pmax is 12 or 14 (see table [TT!), therefore we use a basis set based on 
P — 12 

max 

As well as using the polyad number to truncate the size of the basis, we employ a further 
truncation technique by specifying an upper limit for the eigenvalue calculations, i.e. con- 
struct the basis set such that it provides energy values up to a limit of E'max- This is based 
on an estimation whereby eigenvalues of our ID basis functions are summed together before 
they are considered for matrix element calculations, and the active space is constructed 
using basis functions whose eigenvalues sum together to have E < E'max- For the present 
study we use E^g,^/hc = 10 000 cm~-^. 

These precautions are particularly important for the SO3 molecule, as its larger mass 
(compared to, for example, XH3 systems) gives rise to small rotational constants, which in 
turn requires calculations up to high J value to ensure adequate coverage of transitions for 
a given temperature. This means that any unnecessary basis functions will prove compu- 
tationally expensive. In addition to this basis set minimisation, we can reduce the size of 
the Hamiltonian further by making use of group theory. SO3 has V^i^(M) molecular group 
symmetry and the spin-0 Bosons which make up the constituent atoms of this molecule allow 
the molecular ro- vibrational wavefunctions to exhibit the symmetry of only two of the six 



6 



TABLE II: Convergence of basis set viewed for some vibrational band centres (cm ^) for ^^S-'^^Oa. 







W US. ^ 


p —in 

^ max — -"-^ 


P — 1 9 
^ max — -"-^ 


P —14 
max — -"-^ 
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^1 




1 0^4 Q9 




1065 75 


1 065 74 


















1391 52 
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o,, 0^ 
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1 589 81 


1 587 64 


1 586 46 


1 586 30 


o (J o\ 
07-^4V 4 — ^ / 




1 591 1 


1 587 61 


1 586 43 


1 586 97 


^1 ~1" ^^2 




1 5fiO fiO 


1 565 51 


1 565 '^'^ 


1 565 ^^9 


7/0 -1- Ol/A (1 A 




Loo I .00 


1 ^^f\ '^R 
lOOKJ. 00 


1000. Oc/ 


1 ^f^f^ 47 
1000.^ I 


Z/2 + 2zy4(/4: 


=2) 


1558.52 


1557.12 


1556.45 


1556.37 


2/^2 + T^4{h- 


= 1) 


1525.61 


1524.81 


1524.48 


1524.46 


3z/2 




1492.35 


1449.81 


1490.76 


1490.76 


^2 + ^^3(^3 = 




1884.57 


1881.82 


1881.53 


1881.51 


3^/3(/3 = l) 




4136.39 


4138.88 


4126.78 


4125.92 



irreducible representations of this group in order to satisfy the Pauli Principle; namely the 
A.'^ and A'/ representations. This reduces both the number of Hamiltonian matrices we need 
to consider and, since E symmetry Hamiltonian matrices are larger, their size. 

The two other factors which are important in our spectral calculations are i) the wavenum- 
ber range of the desired synthetic spectrum, and ii) the temperature at which we wish to 
simulate it. The quality of a computed spectrum will become sensitive to E'max as the tem- 
perature increases; we need to ensure that we calculate all energy states that are significantly 
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populated for the given temperature. This can be checked using the temperature-dependent 
partition function: 

Q = Y, giexp{-E,/kT), (10) 

i 

where gi is the total degeneracy of the ro- vibrational state i with energy Ei^ with the sum 
running over all energies at the absolute temperature T, and k is Boltzmann's constant. 
The total degeneracy is given by (2 J + 1) times the nuclear spin degeneracy, which for the 
present case of ^^S^^03 is simply 1 for both the A'^ and A'/ symmetries, given that the nuclear 
spin of ^^O and ^^S are zero. For a given temperature, we can determine the contribution 
of various states to the value of Q. We can then check that Q converges to a specific value 
as Ei tends to infinity; as T increases we require a greater coverage of higher-lying energy 
states. For T = 298.15 K we find that Q converges to better than 1 % at J = 85, with a 
value of Q = 8089.262. Therefore calculations spanning all J's up to 85 should be sufficient 
for simulating spectra at this temperature. Figure [H shows the value of Q as a function of 
all energy levels having J quantum number up to a maximum value, Jmax, at an absolute 
temperature of T = 298.15 K. As we include energies in the summation for increasing values 
of J we see that the associated energy levels contribute less and less to the value of Q, until 
it converges to a limit. 




FIG. 1: Convergence of partition function for different values of Jmax- 



An ab initio PES was previously published by Martin based on coupled cluster meth- 



8 



ods {l^. As an intial test of the TROVE procedure we used the quartic force field by 
Martin to compute fundamental term values, and to test our convergence procedures, using 
a polyad truncation scheme of Pmax= 16. We found discrepancies between our results and 
the values published, particularly in the value of the V2 fundamental term value. We made 
a substitution to the symmetry-adapted force constant F22 associated with this vibration, 
taking a scaled value from a previously published force-field [isl, which resulted in a fair 
improvement. Table [TTT] illustrates the differences between the values and our final converged 
result for the fundamental vibrations. 

TABLE III: A comparison of the fundamental term values (cm~^) between this work and Mar- 
tin Q (computed with TROVE). 



Obs. [8] Martin Martin F22 This work 

z/i 1064.92 1063.36 1064.22 1065.75 

V2 497.57 428.36 487.10 498.48 

z/3 1391.52 1386.81 1386.85 1387.45 

z/4 530.09 527.35 527.32 528.59 



Martin, with a substituted F22 value 



B. Intensity simulations 

The simulation of absorption spectra can be broken down into two main parts: first, 
the eigenfunctions and eigenvalues of the numerically constructed Hamiltonian matrix are 
calculated by a diagonalisation procedure; second, these eigenfunctions are used to compute 
transition dipoles, line strengthes S', Einstein A coefficients and intensities / for allowed 
transitions. For ^^S^^Oa, the rigorous selection rules determining allowed electric dipole 
transitions are A J = J' — J^'' = 0, ±1 [J" + J' > 1), and symmetry selection rules A[ ^ A![. 
The intensity of a transition between given states is given by the formula 

S{f^i). (11) 

where /(/ i) is the transition intensity for a transition from state i with energy Ei to 
state / with energy Ef, with hcvif = Ef — Ei. Q is the partition function. The value of 
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(47reo)3/ic Q(T) 



X 



1 — exp 



if 



kT 



Q must be converged for the absolute temperature T, and with respect to the wavenumber 
range in which we wish to simulate our spectrum. S{f ^ i) is the line strength, which is 
defined by the following integration: 

s{f ^t) = 9nsJ2 E K^L^V^i*^;^)P (12) 

for a transition between initial state i with rovibrational wavefunction and final state 
/ with wavefunction ^iP. Here, ^'ns is the nuclear spin statistical weight factor, and fiA is 
the electronically averaged component of the molecular dipole moment along the space-fixed 
axis A = X^Y^Z. The quantum numbers rrii and nif are projections of the total angular 
momentum J (in units of h) on the laboratory fixed Z axis, for the initial and final states, 
respectively. 



Maki et al js-llll] reported an extensive high-resolution study of a number of fundamental, 
combination and overtone bands of ^^S^^Oa. Their principle aim was to obtain accurate 
wavenumber measurements, but relative intensities were also measured. In the present work 
we convert these data into absolute intensities by normalizing to the theoretical intensities 
obtained with TROVE at T = 298.15 K as described below. 

The measurements available to us cover three spectral regions: 405-708 cm~^ (focusing on 
^^4, 2z/2-^^2, ^^2+^4-^2, ^2+^4-^4, ^1-^4, 2u^^~^'' -u^ and 2iy^^~'^'^ — z/4), 1200-1680 cm~^{iys)^ 
and 2500-3280 cm~-^ (2^/3^^"^^). Each measurement was made at different values of pressure. 
For the 405-708 cm~^ window, the measurements were performed at 0.409 and 2.04 Torr, 
and for the measurements between the 1200-1680 cm~^ window 0.16 Torr and 0.7 Torr 
was used. 560 lines were measured at 0.7 Torr, however 439 of these had relative intensity 
values which were negative. We therefore did not use this higher pressure measurement at 
aU. 

In Table [IV] we compare the numbers of lines identified in each measurement to the 
numbers of lines computed using TROVE. The latter numbers are the subject of the fol- 
lowing selection criteria: J < 85, intensity cut-off, /(/ ^ i) > 10~^^ cm/molecule, and the 
wavenumber window, 0-4000 cm~^. Experimental lines with negative relative intensities 
were left out of the analysis. 

To normalize the experimental intensities, the experimental relative data from each 
spectral window and each different pressure were scaled to match the theoretical values 
computed at T = 298.15 K. The scaling factors obtained through a minimization pro- 
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cedure using all selected experimental lines with non-zero intensity (see Table [Vj) are 
6.571x10-21 cm/molecule (405-708 cm-^ 0.409 Torr), 1.838x10-^1 cm/molecule (405-708 
cm-^ 2.04 Torr), 4.823x10-^° cm/molecule ( 1200-1680 cm-^ 0.16 Torr), and 1.328x10-^1 
cm/molecule (2500-3280 cm--^, 4.99 Torr). With these factors the 'experimental' intensities 
match the theoretical values reasonably good, for example for 405-708 cm-^ the agreement 
is within about 7.02 x 10-^^ cm/molecule and 5.05 x 10-^^ cm/molecule, at 0.409 and 2.04 
Torr, respectively. 

Having the absolute intensities derived, band intensities were estimated for each experi- 
mental band as the sum of individual line intensities. In Table IVl these 'experimental' band 
intensities S^^"^ (cm/molecule) are compared to the theoretical values obtain by summing 
intensities (a) from all TROVE lines from a given window and (b) only from lines with 
experimental counterparts present. This was done separately for each spectral range, and 
each measurement pressure therein. In Table |V] these quantities are referenced to as 5^^^^^^^ 
and S^^^ for the 'total' and 'reduced' band intensities, respectively, and compared to 5^^^^. 
The ratio S^^^ to S^^"^ also shown in Table |V] demonstrates the good quality of the proce- 
dure employed as well as of our dipole moment. For example at 0.409 Torr, the differences 
between 'experimental' and theoretical band intensities are within about 20% for all bands 
from the 405-708 cm--*^ region with the exception of 2iy^^~^'^ - (see also discussion below). 
It should be stressed here that only one scaling factor for all eight bands from this window 
was applied at a given pressure. The difference between two theoretical band intensities 
S^^l^ and S^^^ gives a measure of the missing experimental transition data. According to 
Table |V] even stronger bands miss more that 50% of the total intensity. 

In Table |V] we also show theoretical values of vibrational transition moments defined as 



and ^^^^^ and "if^f^l are the vibrational eigenfunctions of the 'initial' and 'final' states, re- 
spectively, variationally computed using TROVE and /Ha is the component of the molecular 
dipole moment along the molecular- fixed axis a = x^y^ z. 

Figure [2] presents an overview of the simulated spectrum (T = 298.15 K) with TROVE 
and experimental absorption spectra of SO3 for the whole simulation range up to 4000 cm-^. 




(13) 



where 




(14) 
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It reveals the gaps and limitations of the available experimental data. Our intensities based 
on the ab initio DMS are in very good qualitative agreement with experiment. Figure [3] 
shows the 'forbidden' rotational band as a stick spectrum. It should be noted that the 
microwave measurements from Ref. Q| do not have any intensities reported. In Figure [H 
all eight bands from the 405-708 cm~^ region are combined into one graph for each pressure 
to illustrate the quality of the corresponding experimental data. This figure suggests that 
the 0.409 Torr data are generally more reliable. This is also refiected by the ratio values 
S^^^/S^^^ from Table IVl which are significantly closer to 1 at the lower pressure. Of the data 
present for two conditions we therefore place preference on scaled intensity values obtained 
at the lower pressure. Finally, Figure [5] presents a detailed comparison for the all bands 
from the three spectral regions studied in this work in the form of stick diagrams. 



TABLE IV: Comparison of calculated (TROVE) and experimental [8|] band centers and numbers 
of line transitions. 



Band 


Obs. 


Calc. 


Pressure PI" 


Pressure P2'' 


TROVE 




497.57 


498.48 


773 


1265 


5422 


Z/4 - Uq 


530.09 


528.59 


996 


2052 


12195 


Vl - Z/4 


534.83 


537.16 





69 


15147 


Z/2 + - ^2 


530.33 


528.87 


84 


571 


12477 


2V2 - V2 


497.45 


496.88 


112 


704 


7171 


V2 + V/^ -Z^4 


497.81 


498.76 


47 


602 


27182 


iv^ — V\ 


530.36 


528.79 


116 


775 


31096 


ZZ/4 — Z/4 


529.72 


527.91 


39 


455 


13718 


Z/3 - Uq 


1391.52 


1387.45 


2014 




14441 


24^='^ - V, 




2770.29 


1527 




18115 


f - 






25 




3439 



For the z/3 and 2z/3 bands, measurements were taken at 0.16 Torr and 4.99 Torr respectively. 
Pressure values are not recorded for microwave measurements [3]. The remainder are the 
bands within the 405 - 708 cm~^ window, measured at 0.409 Torr (Maki et al 
Bands within the 405-708 cm~^ window, measured at 2.04 Torr (Maki et al 5 
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TABLE V: Vibrational band intensities 5^^^, S^tot^, S^^^ in cm/moleculex 10"-^^, and calculated 
transition moments /2if in Debye. Pi and P2 refer to the different pressure measurements within the 
same wavenumber region (see Table [IV|) . S^^^/S^^"^ is the ratio of the theoretical reduced and total 
band intensities (see text). S^^l^ is the theoretical band intensity computed by summing all TROVE 
lines. S^^^P is the experimental band intensity obtained from a summation over all experimental 
values after scaling factors applied (see text). S^^^ is the theoretical band intensity computed 
using only lines for which experimental counterparts exist. A^red is the number of matched lines. 
Band 5^'^ S^^\^ S^^\^/S'-P N,,^ jl,i/D 







Pi 


P2 




Pi 


P2 


Pi 


P2 


Pi 


P2 








2.987 


1.537 


3.705 


2.915 


1.559 


0.98 


1.01 


773 


1265 


0.158 






4.258 


3.112 


5.949 


4.310 


3.149 


1.01 


1.01 


995 


2052 


0.200 


2V2 - V2 




0.116 


0.470 


0.661 


0.101 


0.411 


0.88 


0.88 


112 


704 


0.221 


1^2 + ^4- 




0.062 


0.322 


0.528 


0.052 


0.251 


0.84 


0.78 


84 


571 


0.199 


l>2 + 1^i- 


V4 


0.026 


0.260 


0.581 


0.022 


0.215 


0.84 


0.83 


47 


602 


0.223 


2^4(14 = 


2) - V4 


0.112 


0.589 


0.873 


0.102 


0.485 


0.91 


0.82 


116 


769 


0.283 


22/4(^4 — 


0) - V4 


0.026 


0.222 


0.405 


0.015 


0.179 


0.57 


0.81 


38 


454 


0.196 


Vi - 1/4 






0.009 


0.101 




0.003 




0.29 




69 


0.039 






39.490 




44.440 


39.490 








2014 




0.321 


2vz 




0.093 




0.119 


0.093 








1527 




0.012 



IV. CONCLUSIONS AND DISCUSSIONS 

The misplaced theoretical bands in Figures [H and [5] indicate that our ab initio PES of 
SO3 requires improvement (see also band centers in Table HVT) . The theoretical U2 frequencies 
have an rms deviation of 0.91 cm~^ when compared with the experimental data, which is 
relatively small compared to the deviation for z/3 of 4.07 cm~^. This is to be expected since 
our PES is purely ab initio computed at a modest level of theory. We are planning to refine 
this surface by fitting to all experimental data available. 

Table |V] outlines the quality of the intensity scaling procedure, in which the relative 
values of the experimental intensities where converted to absolute values (cm/molecule). 
For the 0.409 Torr and 2.04 Torr measurements our comparisons mostly agree to within 
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FIG. 2: Overview of the simulated absorption (T = 298. 15K) spectrum (TROVE) of SO3 compared 
to experiment scaled to the theoretical intensities (see text). 

20%, with the exception of the 2^/4^"^^ - band measured at 0.409 Torr which shows 
nearly a 50% difference, and the Vi - band measured at 2.04 Torr with 80% uncertainty. 
The latter can probably be attributed to both the small number of lines in each case, and 
residual errors in the transition dipole. In the case of the 2vf^~^^ - band the comparison 
with the 2.04 Torr experiment yields a better value for S^fl j S^^^ than for the 0.409 Torr 
measurement, which suggests that the number of lines available at 0.409 Torr is too low (see 
also Figure Sj) The significance of the results presented in Table IVl and illustrated in Figured] 
is that these give an estimation on the quality of the ah initio dipole moment surface as well 
as of the experimental data. Based these numbers we can place a lower estimate on the 
quality uncertainty for our intensities for each band, for example the experiment and theory 
for the V2 and bands agree at least to within 3% for both pressure measurements, while 
it is only 13% for 2v2 - 1^2^ and between 17% -18% for the remaining bands. 

Our complete room-temperature line list for SO3 containing 349 348 513 transitions can 
be accessed online at www. exomol. com in the ExoMol format described in [34] . It includes 
the transition energies, Einstein coefficients A{f ^ i), and absorption intensities estimated 
for T = 298.15 K. Additionally, a list of 10 878 experimental transitions with absolute 
intensities obtained for T = 296 K is included into the supplementary part of this paper in 
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FIG. 3: Rotational absorption band computed for T = 298. 15K, complete up to J=85. 



a form suitable for standard atmospheric and planetary spectroscopic databases. 

Our future work will be focused on the development of a hot line list for ^^S^^Oa for high 
temperature industrial applications as well as for modelling molecular opacity in atmospheres 
of (exo-)planets and cool stars as part of the ExoMol project [35|] (see www.exomol.com). 
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FIG. 4: Comparison plot of TROVE results and the bands of interest measured between 405 - 708 
cm"-*^ by Maki et al , at 0.409 Torr (above) and 2.04 Torr (below). Points are enlarged in some 
cases for clarity. 
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FIG. 5: Fundamental band comparisons between this work (T = 298. 15K, complete up to J=85) 
and Maki et al [sl. The top and middle panels show comparisons for the U2 and z/4 bands, 
respectively, at 0.409 Torr (left) and 2.04 Torr (right). The bottom panel shows a comparison for 
the z/3 band at 0.16 Torr (left) and the 2us band at 4.99 Torr (right). 
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